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The  research  performed  under  this  contract,  during  the  period  12  January 
1983  through  11  January  1984,  can  be  divided  into  two  major  topics;  Finite- 
Difference  synthetic  seismograms  for  SH  waves;  and  array  analysis  of  the  ground 
motions  from  the  1971  San  Fernando  Earthquake. 

In  section  II,  the  accuracy  and  ease  of  application  of  the  finite-difference 
method  for  generating  synthetic  seismograms  of  SH  wave  propagation  in 
cylindrically  symmetric  media  is  discussed.  The  finite-difference  method  has 
the  advantage  that  arbitrary  density  and  velocity  fields  in  the  medium  may  be 
specified.  A  point  source  is  generated  by  a  simple  transformation  of  a  line 
source.  The  accuracy  of  the  finite-difference  seismograms  in  flat  layered  media 
is  confirmed  by  comparison  with  the  Cagniard  De-Hoop  method.  The  finite 
difference  seismograms  also  agree  with  a  previously  untried  dipping-layer  Cagn¬ 
iard  method.  An  earthquake  radiation  pattern  is  approximated  by  introducing  a 
"near  field"  which  has  permanent  displacement  near  the  source. 

In  section  HI,  the  important  new  field  of  strong  ground  motion  modeling  is 
discussed,  where  velocity  profiles  recorded  from  the  1971  San  Fernando  earth¬ 
quake  demonstrate  strong  path  effects  due  to  topography  of  the  basement  sur¬ 
faces.  The  different  waveform  pattern  observed  along  three  long  profiles  implies 
that  valley  or  basin  structures  can  produce  significant  reverberation  and 
amplification  effects.  It  is  suggested  that  the  proper  geometry  of  basins  and 
earthquake  sources  needs  to  be  considered  in  order  to  predict  waveforms  along 
these  wave  paths. 

Velocity  profiles  along  three  local  arrays  in  the  Los  Angeles  basin  suggests 
that,  within  3  km  in  range,  the  waveforms  are  almost  identical.  It  is  also  sug¬ 
gested  that  higher  frequency  waves,  as  represented  by  ground  accelerations, 
are  less  affected  by  these  large  scale  basin  structures. 
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Proflles  of  displacement  traces,  which  are  not  shown  in  this  study,  display 
characteristics  similar  to  those  observed  in  the  velocity  profiles.  That  is,  the 
waveforms  show  good  correlation  with  variations  in  the  basement  surfaces.  How¬ 
ever.  due  to  its  high-frequency  nature,  the  S  phase  is  more  easily  recognized  in 
velocity  records  than  in  displacement  records. 

We  use  the  2D  finite  difference  method  to  compute  synthetic  seismograms 
for  irregular  structures  which  exist  along  the  various  profiles.  Although  the 
structure  used  in  the  calculation  is  rather  simplified  and  the  numerical  code  is 
for  acoustical  media,  the  general  features  of  the  synthetics  agree  quite  well  with 
the  observed  ground  motions.  We  approximate  Heaton's  (1982)  source  model, 
which  has  two  major  areas  of  dislocation  by  point  sources,  at  depths  of  13  km 
and  1  km.  Using  the  relative  timing,  amplitude  and  signal  duration  of  major 
arrivals  observed  in  the  ground  motions  along  both  profile  I  and  profile  II.  we 
suggest  that  the  coherent  S  phase  (the  S  waves)  is  primarily  from  the  first 
source,  at  13  km.  The  shallow  source  (near  1  km)  contributes  little  to  the  sta¬ 
tions  along  the  profile  II.  However,  it  does  produce  strong  surface  waves  for  the 
stations  in  the  San  Fernando  Basin  region  and  also  contributes  to  the  generation 
of  surface  waves  in  the  Los  Angeles  Basin.  This  study  suggests  that  the  observed 
ground  motions  are  indeed  affected  strongly  by  propagation  paths  along 
different  structural  provinces. 


D.  Finite-Difference  Synthetic  Seismograms  for  SH  Waves 

by/  VUaZe,  Don  Helmberger  and  Hob  Clayton 

INTRODUCTION 

The  understanding  of  earthquake  generated  motions  has  improv-  _ 
significantly  in  recent  years.  This  progress  is  due  to  the  refinement  of  methods 
for  generating  synthetic  seismograms  in  the  attempt  to  replicate  observed 
seismograms  and  the  dramatic  improvements  in  computers  as  well  as  an  ever 
growing  collection  of  observations.  The  synthetic  seismogram  incorporates  our 
best  guess  about  the  earthquake  source  properties  and  the  distortions  caused 
by  propagation  through  the  Earth.  Our  knowledge  of  sources  and  earth  struc¬ 
ture  is  dependent  on  the  accuracy  and  versatility  of  the  method  for  synthetic 
seismogram  generation. 

The  earth  was  first  approximated  by  a  wholespace;  it  then  evolved  to  a 
balfspace,  and  it  is  currently  a  stack  of  homogeneous  layers  in  which  the 
number  of  layers  is  controlled  by  one’s  computer  budget.  In  this  paper,  we 
relax  some  of  these  strong  limitations  by  allowing  dipping  structure  and  intro¬ 
ducing  a  procedure  whereby  2-D  finite  difference  (FD)  calculations  can  be 
mapped  into  synthetic  seismograms  with  the  -proper  point  source  shear  disloca¬ 
tion  characteristics.  This  procedure  has  the  advantage  of  allowing  for  arbitrary 
density  and  velocity  fields  in  two  dimensions. 

This  paper  only  treats  SH  waves  (horizontally  polarized  shear  waves)  but  the 
method  may  be  extended  to  P-SV  waves.  We  will  begin  by  reviewing  line  source 
theory,  which  is  the  physical  basis  for  the  FD  method.  The  similarity  of  the  gen¬ 
eralized  ray  theory  (GRT)  expressions  for  strike-slip  and  dip-slip  types  of 


sources  leads  to  identification  of  the  vertical  radiation  patterns  required  for  the 
line-source  FD  source.  FD  seismograms  using  these  radiation  patterns  are  seen 
to  agree  well  with  GRT  seismograms.  Finally,  the  utility  of  FD  methods  with  the 
proposed  mapping  is  shown  for  several  applications. 


SOURCE  REPRESENTATION 

Consider  a  dislocation  source  which  is  situated  below  a  layer  as  shown  Fig¬ 
ure  1.  The  displacement  at  the  free  surface  can  be  obtained  by  summing  gen¬ 
eralized  rays,  after  Helmberger  and  Malone  (1975) 

<« 

where 


v,(t)  =  Vf  5-[ 


(2) 


and  t  is  time  in  seconds,  r  is  horizontal  distance  between  source  and  receiver  in 
km,  n  is  the  number  of  rays  used  to  approximate  the  response,  M0  =  moment  in 
dyne-cm,  p0  =  density  in  grams/cm3  at  source,  D(t)  =  dislocation  history  across 

the  fault  element,  and  -  =  far-fleld  time  function.  A*  and  Ag  contain  the 

dependence  on  azimuth  and  mechanism  and  are 


A*  =  cos  20  cos  X  sin  6  —  4- sin  20  sin  X  sin  26 
Ag  =  -sin  0cos  X  cos  6  -  cos  0sin  X  cos  26 


(3a) 

(3b) 


where  0  =  strike  from  the  end  of  the  fault,  X  =  rake  angle,  and  6  =  dip  angle. 
The  fault  orientation  conventions  are  given  in  Figure  2. 


Receiver 


are  the  strike-slip  and  dip-sup  vertical  radiation  patterns  and  p0  is  the  shear 
wave  velocity  in  km/sec  at  the  source,  e  is  -1  for  receivers  above  the  source, 
but  +1  for  receivers  below  the  source.  is  the  product  of  transmission  and 

reflection  coefficients  for  the  i11*  ray,  p  =  8t*?a  is  the  ray  parameter,  and 

P 

<?.=<;=--  p*>* 

“ i 


*i  -  P  +  Va  ♦  Vi  (  2i  -  l)  h, 


and  (-£),  is  determined  from  tj(p). 
at 

Note  that  V,(t)  and  Vg(t)  do  not  depend  on  azimuth  or  mechanism  but  con¬ 
tain  only  the  vertical  radiation  patterns.  This  separation  allows  the  two  Vj’s  to 
be  constructed  by  equation  (2)  and  the  azimuthal  pattern  and  mechanism  to  be 
calculated  by  equation  (l). 

If  we  allow  fix  =  pg.  equation  (2)  simplifies  for  t*  >>  T,  where  T  is  the  source 
£ 

duration  and  U  =  —  to 

p 


V1  =  2 


•in  a  H(t  -  t,) 


for  a  strike-slip  fault  and 


„  cos  a  H(t  -  t-) 


for  a  dip-slip  fault  where  R  is  the  source  to  receiver  distance  and  a  is  the  angle 
between  the  vertical  and  a  line  connecting  the  source  and  receiver  as  shown  in 
figure  2.  Thus,  the  trigonometric  functions  specify  the  well-known  radiation 
patterns.  Substituting  (5)  into  (l)  and  letting  A  =  0°  and  6  =  90°  for  the  strike- 
slip  case,  we  obtain 


v(r,z.t)  = 


Me 


4trp0 


2F0(  S—")  cos  20  dD^  “  ^  1 


0s 


dt 


7R 


(7) 


where  the  amplitude  of  the  displacement  is  given  in  centimeters,  M0  in  dyne-cm, 
Po  in  grams/cm9  ,  fi  in  km/sec,  R  in  km.  and  F0  =  10_8°  for  unit  conversion  and 
y  =  D(«o)  for  strength  normalization.  We  will  express  most  of  our  results  with 

M0 

- - Fq  =  1  or  for  a  moment  of  Mq  =  4  n  po  x  10*"  dyne-cm. 

4  n  Po 

As  R  and  therefore  t,  becomes  small,  we  need  to  improve  our  approxima¬ 
tions  by  including  a  pseudo-near-field  term,  see  Helmberger  (1983) 


V, 


=  Jr  R®(p)  Im( 


V*  dt7 


which  simplifies  to 


V  2  t  sin  a  H(t  ~  tf) 
l"  fi*  R  R 


and 


v,= 


which  simplifies  to 


v  -  2  t  cos  g  H(t  -  If) 
8  “  p*  R  R 


(Ba) 


(Bb) 


(9a) 


(9b) 


These  are  not  the  complete  near-field  effect!  but  only  compatible  with  the 
asymptotic  nature  of  the  solution.  If  p  were  to  include  the  exact  near-field 
solution  for  the  radiated  SH  waves,  we  would  have  to  include  the  near-field  con¬ 
tributions  of  P  and  SV  as  well,  as  described  in  Helmberger  and  Harkrider  (1978). 

The  solution  for  a  symmetric  line  source  assuming  the  same  geometry  as 
displayed  in  Figure  1  becomes 


* (r.z.t)  =  i-[  f(t)  •  *  ] 


*  =  (ii) 

l«i  Vz 

with  the  same  definitions  given  for  expressions  (l)  and  (2).  A  generalization  of 
line  source  theory  to  include  dipping  structure  is  given  by  Hong  and  Helmberger 


(1977). 


If  we  suppose  pl  -  (3%  as  for  the  point  source  discussed  above,  we  obtain 


2  H(t  -  — ) 


where  the  factor  of  two  is  the  free  surface  effect.  Substituting  (12)  into  (10),  we 
obtain 


$(r,z,t)  = 


S*t- 


which  is  the  well-known  line  source  solution.  By  comparing  expression  (ll)  with 
(2)  and  the  near-field  effects  just  discussed,  we  can  determine  the  constants  and 
mapping  required,  namely  let 
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(14) 

where  $j  is  the  output  of  the  FD  code  with  the  following  source  descriptions 

♦,  =  *„  ( 5%«4 1  -W-V)  -  (15.) 

for  the  strike-slip  orientation  and 

♦,  =  *.(^-)t  (ISb) 

for  the  dip-slip.  The  constant  is 

*o  *  VI  (18> 

The  point  source  seismograms  due  to  shear  dislocations  are  generated  with 
equation  (l)  using  the  Vj  from  equation  (14)  which  use  from  the  FD  code. 

The  separation  of  variables  between  the  Vj‘s  and  the  Aj's  allows  two  different 
ways  to  use  equation  (l).  If  several  synthetic  seismograms  are  desired,  one  cal¬ 
culates  both  Vj's  from  two  FD  runs  and  combines  them  with  the  appropriate  A* 
and  Ag  for  each  seismogram.  If  only  one  seismogram  is  desired,  one  can  com¬ 
bine  the  Vj's  with  the  appropriate  A's  to  create  the  source  and  only  run  the  FD 
code  once. 

Ve  will  next  discuss  the  FD  method  followed  by  comparisons  of  FD  against 
the  line  source  CRT  results  as  given  by  equations  (10)  and  (11),  and  followed  by 
point  source  mapping  comparisons  and  dipping  GRT  results. 


FINITE  DIFFERENCE  TECHNIQUES 


The  finite-difference  (FD)  SH  code  used  is  similar  to  the  acoustic  code 
described  by  Brown  and  Clayton  (1977).  The  code  is  fourth-order  for  regions 
homogeneous  in  velocity  and  second-order  for  regions  heterogeneous  in  velo¬ 
city.  see  Alford  et  al,  1974,  for  definition  and  examples  of  order  for  FD  codes. 
The  fourth-order  scheme  is  used  to  combat  grid  dispersion,  which  slows  down 
the  propagation  of  the  higher  frequency  waves.  Figure  3  shows  the  actual  phase 
velocity  divided  by  the  nominal  phase  velocity  plotted  against  period  in  seconds 
and  wavelength  in  number  of  grid  intervals.  The  line  with  solid  circles  indicates 
the  phase  velocity  of  the  fourth-order  FD  algorithm  used  for  homogeneous 
regions  and  the  line  with  hollow  circles  indicates  the  phase  velocity  of  the 
second-order  algorithm  used  for  heterogeneous  regions.  The  fourth-order 
scheme  allows  propagation  at  approximately  the  correct  velocity  for  energy  with 
wavelengths  down  to  five  grid  points,  or  one-half  the  limit  of  10  points  per 
wavelength  of  a  second-order  scheme.  The  second-order  heterogeneous  condi¬ 
tions  are  used  because  of  the  difficulty  of  deriving  simple,  stable  fourth-order 
heterogeneous  schemes.  The  second-order  heterogeneous  scheme  is  sufficiently 
accurate  for  the  purposes  of  this  paper  and  is  faster. 

Absorbing  boundary  conditions  are  imposed  on  the  sides  and  the  bottom  of 
the  FD  grid  as  described  by  Clayton  and  Engquist  (i960),  and  the  top  of  the  grid 
is  made  a  free  surface  for  SH  waves  (reflection  coefficient  =  1)  by  the  method  of 
Alterman  and  Karal  (1968). 

Several  schemes  for  source  insertion  have  been  described  in  the  literature. 
Imposing  the  displacement  as  a  function  of  time  either  at  one  side  or  at  the  bot¬ 
tom  of  the  box  works  if  the  source  is  too  far  from  the  receiver  to  be  embedded 
in  the  box  (see  Boore,  1972,  for  example).  A  simpler  approach  imposes  an  initial 
displacement  in  a  source  region  for  the  first  two  time  steps  and  leaves  the  grid 
source  free  for  the  remaining  time  steps  (Alford  et  al,  1974).  A  third  method 


commonly  used  is  described  in  Alterman  and  Karal  (1988).  Essentially,  it  solves 
the  wave  propagation  twice  in  the  source  region,  once  imposing  the  source  and 
once  without  a  source,  and  combines  the  results  so  that  the  source  box  does  not 
act  as  a  rigid  reflector.  This  method  has  the  advantage  that  the  source  may  be 
'left  on"  for  as  long  as  desired.  We  use  the  method  of  Alterman  and  Karal 
(1968),  except  that  fourth  instead  of  second  order  boundary  conditions  are  used 
to  match  the  source  region  with  the  surrounding  region.  This  method  is  used 
because  we  desire  a  source  within  the  grid  and  the  wavelengths  of  interest  would 
compel  unacceptably  large  source  boxs  with  the  method  of  Alford  et  al  (1974).  A 
9  by  9  grid  point  source  is  used,  which  is  necessarily  homogeneous  for  our 
source  generation  program. 

Figure  4  illustrates  the  various  sources.  The  frames  in  Figure  4  are 
snapshots  of  displacements  throughout  a  small  grid  at  a  given  time.  The  left 
frames  are  the  displacements  25  timesteps  into  the  calculation  or  when  the 
transient  wave  is  halfway  to  the  edges  and  the  right  frames  are  after  50 
timesteps  or  when  the  transient  wave  has  propagated  almost  to  the  edge  of  the 
box.  The  amplitude  scale  for  each  plot  is  different.  Positive  displacements  are 
filled  in  with  black,  but  negative  displacements  are  difficult  to  discern.  All  three 
types  of  sources  have  a  transient  wave  which  moves  outward,  but  the  sin  6  and 
cos  6  sources  also  have  a  two-lobed  displacement  pattern  which  is  permanent. 
This  permanent  displacement  is  analogous  to  permanent  deformation  caused  by 
movement  on  a  fault. 

The  source  for  Figures  5  through  8  is  symmetric  and  generated  analytically 
by  convolving  an  excitation  function  with  the  impulse  response  for  a  line  source 
which  is  given  in  expression  13.  Figures  9  through  12  use  the  dislocation 
sources  given  in  Figure  4.  These  source  functions  introduce  permanent  dis¬ 
placement  to  the  source  area,  just  as  would  be  expected  from  an  actual 


Figure  4.  Snapshots  of  displacements  for  the  symmetric  source,  the  sin  6 
mechanism  source,  and  the  cos  6  mechanism  source.  The  left  frames 
show  displacements  after  25  timesteps,  the  right  after  50.  Amplitude 
scales  are  not  uniform.  Positive  but  not  negative  displacements  are  filled 
in  with  black. 


LINE  SOURCE 


Figure  5.  Comparison  of  the  line  source  synthetic  seismograms  gen¬ 
erated  by  the  GRT  and  FD  methods.  Results  are  for  a  9  Ion  thick  layer 
over  a  halfspace.  The  source  is  9  km  below  the  layer  and  the  receivers 
are  on  the  surface  and  range  from  40  to  780  km  in  horizontal  distance 
from  the  source.  Both  sets  of  seismograms  are  convolved  with  a  tra¬ 
pezoidal  time  function.  The  amplitudes  are  absolute. 
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Figure  8.  Comparison  of  the  line  source  synthetic  seismograms  gen¬ 
erated  by  the  CRT  and  FD  methods.  Results  are  for  a  9  km  thick  layer 


over  a  halfspace.  Results  are  for  a  layer  which  ranges  in  thickness  from  9 
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km  above  the  source  to  45  km  at  a  range  of  900  km. 
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POINT  SOURCE 
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Figure  7.  Comparison  of  the  point  source  synthetic  seismograms  gen¬ 
erated  by  the  CRT  and  FD  methods  for  the  same  flat-layer  geometry  as  is 
used  for  Figure  S.  Amplitudes  are  absolute  and  may  be  scaled  to 
moment. 


Figure  8.  Comparison  of  the  point  source  synthetic  seismograms  gen¬ 
erated  by  the  GRT  and  FD  methods  for  the  same  dipping-layer  geometry 
as  is  used  for  Figure  8. 
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Figure  9.  Comparison  of  the  point  source  cos  6  mechanism  synthetic 
seismograms  generated  by  the  CRT  and  FD  methods  for  the  same  flat- 
layer  geometry  as  is  used  for  Figure  5. 
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Figure  10.  Comparison  of  the  point  source  sin  6  mechanism  synthetic 
seismograms  generated,  by  the  GRT  and  FD  methods  for  the  same  flat- 
layer  geometry  as  is  used  for  Figure  5. 


DEEP  STRIKE-SLIP  MECHANISM 
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Figure  11.  Comparison  of  the  point  source  cos  6  mechanism  synthetic 
seismograms  generated  by  the  CRT  and  FD  methods.  Results  are  for  the 
same  flat-layer  geometry  as  is  used  for  Figure  5  except  that  the  source  is 
E5  instead  of  9  km  below  the  layer. 


Figure  12.  Comparison  of  the  point  source  sin  6  mechanism  synthetic 
seismograms  generated  by  the  CRT  and  FD  methods.  Results  are  for  the 
same  flat*layer  geometry  as  is  used  for  Figure  5  except  that  the  source  is 
25  instead  of  9  km  below  the  layer. 


dislocation.  For  the  line  source,  the  permanent  displacement  falls  off  with  dis¬ 
tance  as  r_t,  while  the  transient  displacement  falls  off  as  r~1/e. 

Line  source  seismograms  are  described  above,  but  point  source  seismo¬ 
grams  are  desired.  Line  source  seismograms  may  be  transformed  to  point 
source  seismograms  by  equation  14  as  described  above.  A  linear  sum  of  the  two 
sets  of  point  source  seismograms  can  be  used  to  make  seismograms  for  any 
strike,  dip,  and  rake  earthquake  mechanism  for  a  given  velocity  and  density 
structure  as  given  by  equation  1. 

We  use  a  trapezoid  of  3  second  duration  with  1  second  rise  time  and  1 
second  fall  time  as  an  excitation  function  for  the  fiat  layer  models.  A  trapezoid 
of  6  second  duration  with  2  second  rise  and  fall  times  is  used  for  the  dipping 
layer  models.  The  trapezoidal  time  functions  filter  out  frequencies  which  would 
disperse  in  the  grid  and  may  be  thought  of  as  source  time  functions  for  the  syn¬ 
thetic  "shear  explosions"  and  earthquakes  as  is  common  in  modeling  earthquake 
records  (Langston  and  Helmberger,  1975). 

A  symmetric  point  or  line  source  is  easily  grasped  intuitively  if  the  state 
variable  is  pressure  in  an  acoustic  problem,  but  a  symmetric  SH  source  is  more 
difficult.  A  symmetric  SH  line  source  may  be  visualized  as  being  due  to  a  tug  on 
the  line  source.  A  symmetric  SH  point  source,  however,  must  simply  be  con¬ 
sidered  to  have  the  same  SH  amplitude  and  sign  for  every  take-off  angle  and 
azimuth. 

COMPARISON 

A  model  with  one  layer  over  a  halfspace  which  represents  simple  oceanic 
lithosphere  is  chosen  to  compare  the  FD  code  with  the  fiat  layer  GRT  code  which 
has  been  well-tested  (see  Apsel  and  Luco,  1984,  for  example).  The  geometry  is 
shown  in  Figure  2.  The  crust  has  hi  of  9  km  with  a  density  p,  of  2.7  g/cm8  and 


shear  wave  velocity  /?,  of  3.7  km/sec,  and  the  underlying  upper  mantle  has  p2  of 
3.6  g/cms  and  /?2  of  4.8  km/sec.  The  source  is  placed  18  km  below  the  surface, 
or  h2  =  9  km. 

Receivers  in  the  plots  in  Figures  5  through  12  are  positioned  on  the  surface 
at  lateral  distances  of  40  to  760  km  from  the  source.  These  plots  are  reduced  by 
a  4.8  km/sec  velocity. 

Figure  5  contains  synthetic  seismograms  for  the  line  source  problem  with 
symmetric  source.  Results  from  the  GRT  code  are  shown  with  in  Figure  5a  and 
results  from  the  FD  code  are  shown  in  Figure  5b,  where  the  agreement  is  seen  to 
be  excellent. 

Amplitudes  differ  between  the  two  sets  of  synthetics  by  no  more  than  3%. 
and  the  detailed  waveforms  are  nearly  identical.  Even  the  highest  frequencies 
shown  agree  well  for  the  first  portion  of  each  seismogram.  The  later  portion  of 
the  FD  seismograms  show  some  contamination  by  dispersed  high  frequency 
energy.  This  agreement  is  consistent  with  Figure  3.  which  shows  that  energy 
with  periods  above  two  seconds  for  the  homogeneous  regions  and  four  seconds 
for  the  heterogeneous  regions  should  be  propagated  correctly.  The  discrepan¬ 
cies  between  the  traces  are  high  frequency  and  small. 

A  dipping  layer  over  a  halfspace  model  which  represents  a  simple  ocean- 
continent  transition  zone  is  chosen  to  compare  the  FD  code  with  the  dipping- 
layer  GRT  code.  The  crust  is  again  9  km  thick  directly  above  the  source  and 
again  the  source  is  9  km  below  the  crust.  The  thickness  of  the  crust  increases 


, 


linearly  to  45  km  at  a  range  of  900  km.  The  same  geometry  of  receivers  is  used 
as  in  the  previous  case. 

The  results  for  a  symmetric  source  from  the  FD  code  are  shown  in  Figure  6a 


and  from  the  dipping  layer  GRT  code  in  Figure  8b.  The  higher  frequencies  have 
worse  agreement  than  in  the  flat  layer  case,  so  the  excitation  function  used  is  a 
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six  second  trapezoid  with  two  second  rise  and  fall  times.  The  amplitude  of  the 
FD  seismograms  are  3  to  4%  greater  for  the  outer  4  receivers,  but  12  to  18%  less 
for  the  inner  six.  The  discrepancies  between  Figures  8a  and  8b  may  be  due  to 
either  code.  The  most  likely  problems  are  dispersion  in  the  FD  code  as  men¬ 
tioned  above  or  the  necessary  approximation  of  a  dipping  layer  by  a  series  of 
steps  in  the  FD  method.  The  GRT  code  does  not  consider  rays  which  bounce  first 
away  from  the  receivers  before  changing  direction  due  to  the  dip  in  the  layer 
and  travel  to  the  receivers.  The  amplitude  of  these  rays  is  expected  to  be  negli¬ 
gible,  however.  For  the  lower  frequencies  and  the  first  part  of  the  waveform  in 
which  we  are  primarily  interested,  the  agreement  is  excellent. 

The  point  source  synthetics  fare  equally  well.  They  are  generated  by  equa¬ 
tion  14,  described  above.  The  synthetics  are  convolved  with  the  long-period 
WWSSN  instrument  response.  The  longest  and  shortest  period  energy  is 
attenuated  in  transforming  the  line  source  seismograms  to  point  source  seismo¬ 
grams,  primarily  due  to  the  instrument  response.  Figure  7a  shows  traces  gen¬ 
erated  by  the  point-source,  flat-layer  GRT  program  and  Figure  7b  shows  the 
traces  in  Figure  5b  after  the  transformation  described  above.  Similarly,  Figure 
8a  shows  traces  generated  by  the  point-source,  dipping-layer  GRT  program  and 
Figure  8b  shows  the  traces  in  Figure  8b  after  transformation.  The  agreement  in 
waveform  is  excellent  for  the  flat  layer  case  and  good  for  the  dipping  layer  case. 
The  amplitudes  agree  to  within  20%  for  all  cases. 

The  next  question  is  whether  the  FD  seismograms  with  mechanisms 
inserted  agree  with  their  GRT  equivalents.  For  Figures  9  and  10  the  flat  layer 
model  described  above  is  used.  Figure  9  shows  traces  for  the  cos  6  radiation 
pattern,  with  the  GRT-generated  seismograms  in  Figure  9a  and  the  FD  seismo¬ 
grams  in  Figure  9b.  Similarly,  Figure  10  shows  traces  for  the  sin  6  radiation 
pattern.  As  a  further  check,  a  source  25  instead  of  9  km  below  the  layer  is  used 
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in  Figures  11  and  12.  The  agreement  in  waveform  is  excellent,  and  again  the 
amplitudes  agree  to  within  20%.  The  exceptions  are  the  seismograms  recorded 
40  km  away  for  the  source  34  km  deep,  whose  amplitudes  are  substantially  in 
error.  It  appears  that  the  discrepancy  at  A  =  40  km  is  caused  primarily  in 
matching  high  frequency  content  between  the  two  methods.  For  instance,  in 
figure  9  the  GRT  appears  shorter  period  whereas  in  figure  10  the  FD  is  shorter 
period.  These  features  are  probably  caused  by  the  rather  delicate  interference 
between  the  direct  wave  and  first  multiple.  The  interface  at  the  bottom  of  the 
layer  is  actually  smoothed  in  the  DF  code  and  for  angles  of  incidence  less  than 


critical  angle  we  can  expect  a  longer  period  reflection  than  for  GRT  responses. 
Since  we  are  primarily  interested  in  long  range  propagation  in  this  mapping,  we 
can  ignore  the  sharpness  of  the  crust-mantle  transition  which  is  built  into  the 
FD  code  in  this  situation. 

THE  EFFECTS  OF  THE  OCEAN-CONTINENT  TRANSITION  ON  LOVE  WAVES 

If  we  are  to  understand  the  difference  between  the  crust  and  upper  mantle 
under  oceans  and  continents,  it  is  important  to  know  the  geometry  of  the  transi¬ 
tion  zones.  The  sparsity  of  ocean-bottom  seismometers  and  the  relative  quies¬ 
cence  of  passive  margins  precludes  direct  travel-time  inversion  for  velocity 


structure.  Multichannel  seismic  arrays  have  good  resolution,  but  deal  primarily 
with  P-waves  and  can  only  penetrate  the  sedimentary  column  and  occasionally 
locate  the  Moho  discontinuity  (e.g.  Talwani  et  al,  1982,  and  Avedik  et  al,  1982). 
Surface  wave  dispersion  is  often  used  to  examine  velocity  structure,  but  most 
methods  of  surface  wave  dispersion  analysis  require  either  flat-  or  dipping¬ 
layered  models.  Velocity  structure  is  known  not  to  follow  either  flat-  or  dipping- 
layered  geometries  near  ocean-continent  boundaries. 
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The  FD  method  allows  more  general  models  to  be  incorporated  in  synthetic 
seismogram  generation,  as  described  by  Drake  (1972),  Drake  and  Bolt  (1980) 
and  Martel  (1980).  Such  studies  are  two-dimensional,  however,  and  are  more 
complicated  than  the  FD  method  described  in  this  paper.  Some  three- 
dimensional  finite-element  models  have  been  investigated  by  Gregersen  and 
Alsop  (1974)  and  Schlue  (1981),  but  these  models  are  not  very  realistic. 

Several  recent  papers  use  approaches  similar  in  some  ways  to  ours.  Fuyuki 
and  Matsumoto  (1980)  use  a  two-dimensional  P-SV  FD  scheme  to  investigate  the 
effect  of  a  trench  on  Rayleigh  wave  propagation.  Stevens  (1983)  investigates  the 
effect  on  seismic  sections  of  bump  on  the  seafloor  with  a  three-dimensional  P-SV 
FD  code  which  assumes  cylindrical  symmetry.  Kelly  (1983)  investigates  the 
effect  of  near-surface  geological  features  on  the  propagation  of  Love  waves  with 
a  two-dimensional  SH  FD  code. 

In  this  section,  the  effect  of  differing  ocean-continent  transition  geometry 
on  Love  wave  waveforms  is  examined.  Figure  13  shows  the  three  geometries 
used  to  generate  synthetic  seismograms.  In  all  cases,  the  upper  layer  has  a 
density  of  2.7  g/cm8  and  shear  wave  velocity  of  3.7  km/sec  and  the  lower  layer 
has  a  density  of  3.6  g/cm8  and  shear  wave  velocity  of  4.8  km/sec.  Symmetric 
shear  sources  are  used;  these  create  seismograms  which  are  indistinguishable 
from  strike-slip  mechanisms  for  these  ranges.  The  seismograms  for  the 
geometries  shown  in  Figures  13a,  13b,  and  13c  are  presented  in  Figures  14a, 
14b,  and  14c,  respectively. 

The  seismograms  for  different  geometries  are  distinguishable  by  waveform 
and  timing  but  not  by  amplitude  in  synthetic  seismograms.  The  differences  in 
amplitude  are  of  a  similar  magnitude  to  the  errors  noted  above  and  are  not 
clearly  systematic.  The  waveforms  at  the  first  4  receivers,  6  =  40,  120,  200,  and 
280  km  are  identical  because  the  structure  between  the  source  and  receiver  is 
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Figure  13.  Shown  are  the  three  geometries  which  are  used  to  generate 
the  synthetic  seismograms  in  Figures  14a,  14b  and  14c.  The  stars  indi¬ 
cate  sources  and  the  triangles  indicate  receivers.  The  crust  has  fl  =  3.7 
km/sec  and  p  =  2.7  g/cc  and  the  uppermost  mantle  has  0  *  4.8  km/sec 
and  p  *  3.2  g/cc.  The  crustal  thickness  in  each  model  varies  from  5  km 
on  the  left  to  45  km  on  the  right. 


OCEAN-CONTINENT  TRANSITION 


Figure  14.  Suites  of  synthetic  seismograms  are  shown  for  the  3 
geometries  depicted  in  Figure  13.  The  seismograms  correspond  to  dis¬ 
tances  which  range  from  40  to  760  km.  Amplitudes  are  absolute  and  may 
be  converted  to  displacements  when  an  earthquake  moment  is  specified. 


identical.  The  seismograms  recorded  directly  above  the  transition  zone  at  6  = 
360,  440,  and  520  km  are  very  similar.  The  furthest  seismograms  recorded  at  6 
—  800,  680,  and  760  km  are  the  most  distinguishable.  This  pattern  is  fortunate 
because  most  of  the  available  data  is,  of  course,  recorded  on  the  continents. 
The  energy  arrives  slightly  later  in  the  wide  ramp  model  than  in  the  narrow 
ramp  model  which  in  turn  is  slightly  slower  than  in  the  step  model.  This  trend 
may  be  due  to  energy  entering  the  crust  at  points  A  in  Figure  13.  The  slow  cru¬ 
stal  velocity  then  results  in  the  differential  delay  mentioned  above. 

It  has  been  noticed  that  observations  are  more  long-period  than  synthetic 
seismograms  constructed  with  impulsive  sources.  This  effect  may  be  due  to 
drawn-out  smooth  earthquake  source  time  functions  or  it  may  be  due  to  Q 
(attenuation),  which  damps  higher  frequencies  preferentially.  This  effect  means 
that  the  longer  period  shape  of  the  sections  are  the  most  helpful  in  the  modeling 
of  actual  seismograms.  Grand  and  Helmberger  (1983)  find  empirically  that  the 
average  crustal  thickness  between  source  and  receiver  seems  to  determine  the 
shape  of  the  long-period  front  of  the  seismogram.  If  this  idea  works,  then  the 
front  of  the  seismograms  at  6  =  6B0  and  760  km  should  be  similar.  They  are  not 
identical,  but  the  leading  portion  of  the  0  =  760  km  seismograms  for  the  two 
ramp  geometries  are  very  similar.  If  this  idea  proves  usable,  then  one  could 
determine  the  average  thickness  of  the  crust  from  the  front  of  the  seismogram 
and  more  detailed  information  about  the  crust-mantle  transition  from  the  later 
portion  of  the  seismogram,  which  is  perhaps  more  difficult  to  model. 
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III.  Array  Analysis  of  the  Ground  Motions  from  the  1971 
San  Fernando  California  Earthquake 

by  Hsui-Lin  Liu 
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Introduction 

An  important  new  field  of  strong  ground  motion  research  is  modeling  of  an 
earthquake  source  and  the  associated  wave  propagation  in  a  deterministic  way. 
This  type  of  modeling  has  significantly  increased  our  understanding  of  the 
longer-period  parts  of  strong  ground  motions.  However,  a  proper  model  is  often 
difficult  to  obtain.  An  accurate  characterization  of  the  wave  phenomenon.  (  e.g., 
wave  types,  coherent  phase  velocities,  propagation  direction  and  the  wave  varia¬ 
bility  )  is  essential  to  evaluate  the  validity  of  the  models. 

Hanks  (  1975  )  analyzed  the  ground  displacements  from  the  1971  San  Fer¬ 
nando  strong-motion  data  and  found  strong  coherence  appeared  in  the  longer- 
period  waves.  He  also  identified  direct  shear  phases  and  subsequent  surface 
waves  at  certain  ranges.  However,  a  deterministic  source  model  for  the  1971 
San  Fernando  event  is  still  uncertain  as  pointed  out  by  Heaton  (  1982  ).  Heaton 
constructed  models  consisting  of  finite  faults  within  a  half-space  to  explain 
waveforms  observed  at  5  strong-motion  stations  as  well  as  long-period  telese- 
ismic  body  waveforms.  His  preferred  finite  source  model  implies  a  fairly  com¬ 
plex  double  source  earthquake,  with  two  events  of  approximately  equal  size  and 
separated  by  about  A  seconds.  However,  modeling  of  this  type  is  limited  by  our 
incomplete  knowledge  of  the  strong  motions  propagating  through  the  complex 


structure  encountered  in  the  real  earth.  Thus,  we  believe  that  there  is  a  great 
need  for  a  detailed  description  of  the  wave  types  represented  by  the  ground 
motions. 

In  this  study,  we  extend  Hanks’  study  (  1975  )  to  examine  the  ground  veloci¬ 
ties  and  accelerations  observed  by  the  accelerograph  array  during  the  1971  San 
Fernando  earthquake.  A  display  of  the  ground  velocity  profiles  together  with  the 
topographic  features  and  geologic  structure  (  Yerkes  et  al.,  1965  )  will  be 
presented  and  possible  wave  types  and  apparent  phase  velocities  will  be  dis¬ 
cussed. 

Given  the  distribution  of  strong  motion  recordings  of  the  San  Fernando 
earthquake,  it  is  possible  to  construct  profiles  along  three  azimuths.  Figure  3-1 
shows  the  station  distribution  and  corresponding  profiles.  All  profiles  originate 
in  the  epicentral  region.  Profile  I  extends  85  km  southward  across  the  San  Fer¬ 
nando  and  Los  Angeles  basins  to  a  station  on  the  Palos  Verdes  Peninsula.  Profile 
II  extends  95  km  S40*E  along  the  front  of  the  San  Gabriel  mountains  and  then 
across  the  San  Gabriel  and  Los  Angeles  basins.  Profile  111  runs  90  km  N40*W 
across  the  San  Gabriel  and  Tehachapi  mountains. 

The  San  Fernando  earthquake  was  also  well  recorded  by  strong  motion 
instruments  which  are  required  in  high-rise  buildings.  These  buildings  are 
clustered  in  three  locations  on  the  north  side  of  the  Los  Angeles  basin;  the  Mira¬ 
cle  Mile  area  of  Wilshire  Boulevard  {  local  array  1  ),  downtown  Los  Angeles  ( local 
array  2  ),  and  Century  City  (  local  array  3  ).  Station  locations  for  these  dense 
local  arrays  are  shown  in  Figure  3-2. 

Time  histories  of  ground  velocity  and  acceleration  are  ail  taken  directly 
from  EERL,  Caltech  reports  (  1974  ).  Ground  motions  are  rotated  into  radial, 


-35- 


^EriCENTeA 

0(114 

srres  not  in  stuov  i 

1  0(110  1 

ACCEL EROORARN 

MAJOR FAULTS 

<Z%P 

EXPOSED  BASEMENT  ROCK 


6  - n  — 

*  BASEMENT  ELEVATION  IN 
THOUSANOS  or  FEET 
SEA  LEVEL  DATUM 


V  \ 

SANTA  JUSAMMA  MTS 


SANTA  MON/CA  MTS 


mts  \ 


©AREA  I  ©AREA 2  ©AREAS 


1  075 

COSO 

lost 

EM3 

roe 

INI 

toe 

f  ON 

I  134 

R2I7 

C  112 

I 131 

SS5 

1 1ST 

■  IN 

S2K 

1  ISO 

J  Id 

j  IS} 

PACIFIC 

“l 


V.*s 


-r 


OCEAN 


Figure  3*1  The  long  profiles,  I.  11  and  Ill,  and  site  locations  of  accelerometer 
recordings  of  the  1971  San  Fernando  earthquake  (  modified  from 
Hanks,  1975  ). 


transverse,  and  vertical  components  and  are  then  displayed  as  functions  of  time 
and  distance  from  the  epicenter  as  reported  by  Allen  et  al.  (  1973  ).  Since  abso¬ 
lute  time  is  not  available  for  any  records,  some  assumptions  need  to  be  made  to 
correlate  phases  from  one  station  to  another.  We  shifted  records  such  that  the 
apparent  first  shear-wave  arrivals,  which  are  named  S!  by  Hanks  (  1975  ),  are 
aligned  vertically  in  the  profile.  S!  minus  trigger  times  are  listed  in  Table  3-1  for 
all  stations  used  in  this  study.  Unfortunately,  in  some  cases,  the  identification 
of  S!  is  very  difficult,  and  thus  some  of  the  more  distant  stations  are  probably 
not  aligned  on  the  first  shear-wave  arrival.  Once  the  profiles  are  constructed,  we 
find  that  the  ground  velocity  waveforms  are  surprisingly  coherent  from  station 
to  staion,  allowing  us  to  estimate  apparent  phase  velocities,  to  identify  possible 
wave  types,  and  qualitatively  separate  the  source  from  the  path  effects.  The 
features  of  each  profile  will  be  described  in  more  detail  as  follows. 

Long  Profiles 

a.  Profile  I.  The  velocity  traces  together  with  the  corresponding  free  sur¬ 
face  and  basement  surface  topography  profiles  are  displayed  in  Figure  3-3. 
Velocity  traces  begin  with  the  S!  arrivals.  The  two-fault  system  suggested  by 
Heaton  (1982)  is  also  illustrated  in  the  same  Figure.  Hi  is  the  first  hypocenter 
which  ruptured  along  the  Sierra  Madre  fault,  and  4  seconds  later,  the  second 
hypocenter  H2  initiated  on  the  San  Fernando  fault.  According  to  Heaton's 
(1982)  interpretation,  the  high  peak  velocity  observed  at  station  C041  (Pacoima 
Dam)  is  caused  by  rupture  directivity  along  the  Sierra  Madre  fault.  Although  the 
faulting  process  is  rather  complicated,  the  total  source  time  duration  is  only 
about  7  seconds.  We  believe  that  the  signed  duration  at  Pacoima  Dam  (C041) 
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PROFILE  I 


Figure  3-3  Transverse,  radial  and  vertical  components  of  ground  velocity  along 
profile  I.  The  corresponding  free  surface  and  basement  surface  are 
shown  on  the  top.  Dashed  lines  crossing  30  km  to  40  km  indicate  the 
possible  phase  arrivals  of  the  surface  waves. 
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represents  the  approximate  duration  of  the  source. 

One  of  the  most  impressive  features  seen  in  Figure  3-3  is  the  correlation  of 
ground  motion  waveform  with  the  topography  of  the  subsurface  basins.  There 
are  three  stations  C04B,  J145  and  KH5  located  within  the  San  Fernando  Valley. 
The  waveforms  recorded  at  these  stations  are  poorly  correlated  and  the  signal 
durations  are  about  30  seconds.  These  durations  are  much  longer  than  those 
seen  at  neighboring  and  more  distant  stations  located  near  Santa  Monica  Mts., 
such  as  L166,  D068  and  D0S7.  This  contrast  can  easily  be  seen  in  Figure  3-3. 
Coherent  waveforms  begin  to  appear  at  station  L138  and  signal  durations  drop  to 
only  about  10  seconds.  Furthermore,  the  amplitudes  also  drop  quite  drastically 
by  about  one-half  after  passing  the  San  Fernando  Valley  stations.  Amplitudes 
again  increase  at  stations  D057  and  S262  when  the  waves  reach  the  Los  Angeles 
Basin.  The  waveforms  are  quite  coherent  in  the  beginning  6  seconds  and  there 
appears  to  be  little  moveout  from  the  S!  phase.  The  dashed  lines  passing  sta¬ 
tions  L166.  D068,  D057  and  S262  indicate  the  possible  arrival  of  surface  waves, 
which  are  characterized  by  gradual  moveout  from  the  S  phases.  They  also 
appear  dispersive  in  nature,  this  was  also  suggested  by  Hanks  (  1975  )  from  bis 
study  of  particle  trajectories.  It  seems  that,  at  epicentral  distances  of  near  30 
km.  body  waves  dominate  the  waveforms  at  the  bedrock  sites,  and  the  surface 
waves  start  to  develop  rapidly  from  there.  For  epicentral  ranges  beyond  40  km. 
surface  waves  start  to  dominate  the  wavetrain.  At  ranges  of  50  km  and  beyond, 
it  is  hard  to  identify  the  body  waves.  Most  of  the  signals  appear  to  be  surface 
waves  with  visibly  dispersive  character.  The  apparent  moveout  velocity  of  the 
surface  wave  phase  from  S!  phase  is  estimated  to  be  5.7  km/sec.  If  we  assume 
that  the  phase  velocity  corresponding  to  the  S!  phase  is  3.5  km/sec,  then  the 
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phase  velocity  of  these  surface  waves  is  about  2.2  km/sec. 

It  is  useful  to  contrast  the  nature  of  wave  propagation  for  waves  observed  in 
the  San  Fernando  valley  and  the  Los  Angeles  basin.  Peak  ground  velocities 
observed  in  both  basins  are  high  with  respect  to  sites  in  the  Santa  Monica  moun¬ 
tains.  Although  long-period  surface  waves  are  clearly  visible  in  both  basins,  the 
duration  of  these  surface  waves  is  up  to  20  seconds  longer  in  the  San  Fernando 
valley,  the  basin  which  is  closer  to  the  earthquake  source.  The  earthquake  rup¬ 
tured  into  the  northern  part  of  the  San  Fernando  valley  which  is  a  closed  sedi¬ 
mentary  basin.  Thus  it  is  likely  that  the  long  signal  durations  are  caused  by  sur¬ 
face  waves  which  are  trapped  within  the  valley,  thereby  causing  severe  ringing. 
The  short  signal  durations  seen  at  adjacent  stations  indicates  that  these  rever¬ 
berating  waves  which  developed  in  the  San  Fernando  valley  did  not  propagate 
away  from  the  valley.  The  interpretation  from  the  numerical  simulation  result 
will  be  presented  later. 

In  Figure  3-4,  we  show  the  radial  components  of  the  ground  acceleration  for 
stations  located  along  profile  I.  Traces  begin  at  their  trigger  time  and  are 
aligned  vertically  with  the  SI  phase.  Although  there  is  considerable  variation  in 
the  observed  waveforms  from  station  to  station,  the  duration  of  the  high- 
frequency  motions  is  a  fairly  constant  10  seconds.  Furthermore,  strong  phases 
can  be  seen  on  many  records  at  about  4  seconds  and  8  seconds  after  the  initial  S 
arrivals.  We  speculate  that  the  overall  duration  and  timing  of  arrivals  seen  on 
the  relatively  high-frequency  acceleration  time  history  is  directly  related  to  the 
details  of  the  faulting  process.  It  is  suggested  that  the  basin  geometry  has  less 
effect  on  these  higher  frequency  waves. 


along  profile  I. 


Figure  3-4  The  radial  component  of  ground  acceleration 


b.  Profile  II.  Profile  II  is  roughly  along  the  southwest  border  of  the  San 
Gabriel  mountains  and  part  of  the  Los  Angeles  basin.  The  velocity  traces  (begin¬ 
ning  at  S!)  together  with  the  path  profile  are  shown  in  Figure  3-5.  Somewhat 
different  features  are  observed  along  this  profile.  The  signal  duration  of  what 
appear  to  be  mainly  the  SI  phase  is  about  B  seconds  at  station  G110.  Surface 
waves  develop  gradually  with  distance  from  there.  The  amplitudes  along  this 
path  are,  on  average,  lower  than  those  in  Profile  1.  Perhaps  part  of  the  reason 
for  this  is  the  lack  of  thick  sediments  along  profile  II  from  about  30  km  to  50  km. 
Well  dispersed  surface  waves  can  hardly  be  identified  at  this  portion  of  the 
profile.  Instead,  the  S!  phase,  which  is  characterized  by  zero  moveout  can  still 
easily  be  recognized.  Even  at  station  N186,  which  has  an  epicentral  distance  of 
52  km.  the  high-frequency  S!  phase  can  still  be  identified,  this  suggests  that  the 
relatively  thin  sediments  along  profile  II  prohibit  the  rapid  development  of  the 
surface  waves  seen  in  profile  1.  However,  there  may  still  be  a  strong  fundamen¬ 
tal  Rayleigh  wave  which  is  hard  to  distinguish  from  the  S  wave  (  see  Heaton  and 
Helmberger,  1979  ).  The  dashed  lines  crossing  station  G110  to  station  N1B6  indi¬ 
cate  the  possible  surface  wave  phases.  If  we  assume  a  horizontal  phase  velocity 
of  3.5  km/sec  for  the  SI  phase,  then  the  phase  velocity  corresponding  to  the 
dashed  line  is  approximately  2.7  km/sec. 

In  Figure  3-6,  we  show  the  radial  components  of  the  ground  acceleration  for 
stations  located  along  Profile  II.  Peak  amplitudes  and  the  duration  of  high- 
frequency  shaking  are  remarkably  similar  to  that  seen  along  Profile  I.  However, 
the  strong  phases  seen  at  4  and  8  seconds  in  S  wavetrain  on  Profile  I  are  not  as 
evident  in  Profile  II.  Instead,  the  strongest  phase  seen  along  Profile  n  appears  to 
be  about  6  seconds  after  the  beginning  of  the  S  wavetrain. 
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Figure  3-8  The  radial  component  of  ground  acceleration  along  profile  11. 
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e.  Profile  m.  Profile  III  runs  N40°  W  across  mostly  mountainous  terrain  and 
ground  velocity  traces  along  this  profile  are  shown  in  Figure  3-7.  Ground  veloci¬ 
ties  are  strikingly  different  along  this  profile  from  the  other  two  profiles,  which 
cross  deep  sedimentary  basins.  The  duration  and  amplitude  of  the  ground  velo¬ 
city  seem  to  be  considerably  less  on  this  northern  profile.  Heaton  and  Helm- 
berger  (  1979  )  shows  that  displacements  observed  at  station  J142  can  be  rea¬ 
sonably  modeled  with  a  simple  half-space  structure.  This  is  obviously  not  the 
case  for  many  of  the  stations  located  along  profiles  I  and  II. 

Radial  ground  accelerations  observed  along  Profile  III  are  shown  in  Figure 
3-8.  Ground  accelerations  appear  to  be  shorter-period  and  of  shorter  duration 
than  those  observed  along  the  other  two  profiles.  Furthermore,  it  is  difficult  to 
identify  any  coherent  phases. 

In  Figure  3-9,  we  compare  peak  radial  velocities  amplitudes  along  the  three 
profiles.  It  seems  clear  that  peak  velocities  are  larger  on  stations  to  the  south, 
except  for  those  stations  on  or  adjacent  to  the  Santa  Monica  or  San  Gabriel 
mountains.  It  appears  likely  that  the  deep  sedimentary  basins  are  the  principal 
reason  for  this  amplitude  variation. 

In  Figure  3-10,  we  compare  peak  radial  acceleration  amplitudes  along  the 
profiles.  Peak  accelerations  are  surprisingly  similar  for  these  three  profiles. 
This  conclusion  was  previously  reported  by  Trifunic  (  1976  ).  Using  data  mainly 
from  the  San  Fernando  earthquake,  he  concluded  that  peak  ground  motions 
were  greater  on  'soft'  sites  than  on  'hard*  sites,  with  the  effect  being  consider¬ 
able  at  long  periods  and  almost  negligible  at  high-frequencies. 


Log  Velocity  (R)  (cm /sec) 


Figure  3-9  Peak  radial  velocities  vs.  epicentral  distance  for  the  three  long 
profiles. 
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Figure  3-10  Peak  radial  accelerations  vs.  epicentral  distance  for  three  long 
profiles,  1,  II,  and  III. 


Local  Arrays 

In  Figure  3-11,  3-12  and  3-13,  we  show  three  components  of  ground  velocity 
as  funotions  of  time  and  epicentral  distance  for  local  arrays  1,  2  and  3,  respec¬ 
tively.  Since  stations  in  these  arrays  are  not  generally  aligned  along  a  radial 
line  from  the  epicenter,  the  distances  shown  in  Figure  3-11  through  3-13  may 
not  be  indicative  of  the  true  inter-station  spacing. 

The  most  obvious  feature  of  the  ground  velocity  is  the  remarkable  degree  of 
coherence  across  each  array.  Since  the  velocity  waveforms  are  dominated  by 
waves  with  frequencies  of  less  than  5  Hz,  we  expect  that  most  of  the  motion  has 
horizontal  wavelengths  of  greater  than  0.5  km.  Thus,  there  is  good  reason  to 
expect  the  ground  velocity  to  be  coherent  over  these  small  arrays.  However, 
there  are  several  other  implications  of  this  coherence.  First,  the  record  process¬ 
ing  must  introduce  very  little  noise  into  the  velocity  traces.  This  is  especially 
impressive  when  one  considers  that  the  horizontal  components  require  com¬ 
ponent  rotations  before  they  can  be  compared.  Second,  the  effect  of  building 
interaction  is  apparently  minimal  on  velocity  waveforms.  The  average  standard 
derivation  of  peak  velocity  within  each  array  is  about  12%.  Furthermore,  it  is 
difficult  to  see  the  systematic  dependence  of  peak  velocity  with  building  height 
(see  Table  3-1)  within  each  array  that  is  suggested  by  Boore  et  al.  (  1980  ). 

We  show  the  radial  component  of  ground  acceleration  for  local  array  1  in 
Figure  3-14.  Although  less  coherent  than  the  velocities,  accelerations  still  show 
strong  coherence  along  the  profile.  Prominent  phases  can  be  seen  at  2.  4.  and  B 
seconds  after  the  initial  S!  phase.  It  is  difficult  to  see  moveout  of  any  phase  for 
these  local  arrays.  This  is  not  surprising  since  a  variation  in  phase  velocity 
between  4  and  2  km/sec  yields  a  moveout  of  only  .25  sec/km.  Such  small  time 


Figure  3-12  Transverse,  radial  and  vertical  components  of  ground  velocity  in  lo¬ 
cal  array-2. 
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Figure  3-14  The  radial  component  of  ground  acceleration  in  local  array- 
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shifts  could  only  be  detected  through  very  careful  analysis. 

Simulation  by  Finite  Difference  Method 

The  strong  correlation  of  observed  ground  motions  with  geological  struc¬ 
ture  profiles  encourages  our  further  investigation  of  these  irregular  structure 
effects.  We  use  the  acoustical  2D  finite  difference  method  (  Grown  and  Clayton, 
1983)  to  generate  synthetic  seismograms  along  both  profile  I  and  profile  II.  The 
main  emphasis  here  will  be  on  phenomena  which  occur  on  waves  propagating 
through  these  irregular  structures  from  different  source  depths.  We  investigate 
amplitude  amplification,  signal  duration  and  wave  type  development. 

According  to  Heaton's  (  1982  )  double  event  solution  for  the  1971  San  Fer¬ 
nando  earthquake,  there  are  mainly  two  areas  in  which  the  dislocation  is  con¬ 
centrated.  One  area  is  near  the  first  hypocenter  at  a  depth  of  about  13  km  and 
has  a  maximum  dislocation  of  about  2.5  meters.  Another  is  very  near  the  surface 
(  at  depth  of  about  1  km  )  along  the  San  Fernando  fault  with  a  maximum  disloca¬ 
tion  of  about  3.5  meters,  as  illustrated  in  Figure  3-15.  In  order  to  study  the 
effect  of  sources  at  different  depths,  we  calculated  the  structure  response  along 
each  profile  for  three  source  depths,  namely  1  km  .  7km  and  13  km.  The 
difference  method  is  applied  to  the  2D  acoustical  wave  equation 

,  where  z.z  are  the  2D  coordinates  and  v{x,z)  is  the  medium  velocity.  The  nor¬ 
mal  stress  or  the  pressure  P(x,z)  is  continuous  along  all  boundaries  and  van¬ 
ishes  on  the  free  surface.  The  details  are  described  in  a  paper  by  Brown  and 


Figure  3-15  The  dislocation  distribution  along  the  two  faults  in  1971  San  Fernan¬ 
do  earthquake  determined  by  the  waveform  modeling  of  Heaton 
(  after  Heaton,  1982  ). 


Clayton  (  1983  ).  We  used  a  0.2  km  grid  size  for  both  horizontal  and  vertical 
ranges,  this  corresponds  to  about  10  samples  per  wavelength  at  a  period  of  1 
second.  Since  the  maximum  velocity  in  our  model  is  3.5  km/sec,  we  chose  the 
difference  time  step  dt  to  be  0.03  seconds  to  satisfy  the  stability  criteria. 


where  dh  is  the  grid  size.  In  this  way  we  constrain  the  time  step  to  be  small 
enough  that  disturbances  do  not  propagate  faster  than  the  maximum  velocity 
within  a  grid. 

a.  Irregular  Structure  along  Profile  I  We  approximate  earth  structure  along 
our  profiles  as  a  simple  low-velocity  (2.1  Icm/sec  )  basin  embedded  within  a 
high-velocity  (  3.5  km/sec  )  bedrock  half-space,  as  illustrated  in  the  left  of  Fig¬ 
ure  3-16(a).  The  source  is  introduced  as  an  initial  condition  whereby  wave 
motions  in  the  source  region  are  specified  to  be  the  derivative  of  a  Gaussian 
function.  The  details  can  be  found  in  Brown  and  Clayton  (  1983  ).  All  the  profiles 
are  computed  assuming  a  400x300  grid,  which  corresponds  to  dimensions  of 
BOfcmxGOfcm.  Synthetic  seismograms  computed  along  profiles  50  km  long  are 
shown  in  Figures  3-18(a),  3-16(b)  and  3-18(c)  for  source  depths  of  1  km  (  S3  ).  7 
km  (  S2  )  and  13  km  (  Si  ),  respectively.  The  total  signal  duration  is  45-second 
for  all  the  synthetic  traces  presented  in  this  Chapter.  Full  scale  amplitude  is 
0.07  for  all  the  profile  traces.  Thus,  waveform  amplitudes  are  comparable  for 
different  sources. 

Although  the  synthetics  are  computed  assuming  a  2D  acoustic  line  source, 
the  resulting  synthetics  show  good  resemblance  to  the  observed  ground 
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velocttiea  along  profile  1.  In  our  model,  the  San  Fernando  Beam  la  an  obstacle, 
creating  a  shadow  for  those  receiver  sites  on  the  Santa  Monica  Mountains.  The 
strong  amplification  in  the  San  Fernando  Basin  can  be  seen  very  well  for  all 
three  source  depths.  Strong  energy  focusing  is  seen  near  the  edge  of  the  basin. 
This  phenomenon  is  discussed  in  more  detail  in  Appendix  II  and  is  illustrated  for 
a  simpler  geometry  in  Figures  11-5,  II-6  and  II-7. 

Little  energy  from  the  shallow  source  transmits  into  the  ridge  of  basement 
(Santa  Monica  Mountains),  which  separates  the  two  low-velocity  basins.  How¬ 
ever,  for  deeper  source  (  say  at  13  km  )  direct  arrivals  actually  become  stronger 
in  the  ridge  of  basement  since  waves  can  reach  this  region  from  below  the  basin. 

Entering  the  Los  Angeles  Basin,  we  see  that  surface  waves  become  the  dom¬ 
inant  arrivals.  Wave  development  here  is  similar  to  examples  of  down-dip  wave 
propagation  discussed  in  Appendix  n  (  see  Figures  11-1,  II-2,  and  II-3  ).  For  the 
deeper  sources,  energy  reaches  the  Los  Angeles  Basin  without  strong  effects 
from  the  San  Fernando  basin.  This  can  be  recognized  by  comparing  between  the 
profiles  in  Figure  3-16(a)  and  Figure  II-3.  However,  in  this  range,  the  interpreta¬ 
tion  of  wave  propagating  through  the  ridge  from  the  shallow  source  is  not 
straightforward.  We  demonstrate  this  by  comparing  the  synthetics  along  two 
different  ridge  structures.  In  Figure  3-17(a),  the  San  Fernando  basin  is  simu¬ 
lated  by  a  simple  soft  layer  with  an  up-dip  thinning  edge.  The  synthetic  seismo¬ 
grams  suggest  that  little  energy  propagates  through  the  basin  and  into  the  half¬ 
space  region  beyond.  However,  when  a  second  basin  is  added  to  form  a  ridge  as 
illustrated  in  Figure  3-17(b),  then  surface  waves  gradually  develop  in  the  second 
basin  after  passing  the  ridge  area.  Thus,  the  energy  leaking  through  the  up-dip 
edge  of  the  San  Fernando  basin  is  trapped  again  by  the  down-dip  slope  in  the  Los 
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Figure  3-17(a)  Ridge  structure  simulation;  up-dip  edge  of  the  basin.  This  illus¬ 
tration  shows  that  little  energy  can  pass  the  up-dip  basin  edge  and 
reach  the  surface  receivers.  Energy  that  leaks  through  the  slope  pro¬ 
pagates  mainly  toward  the  half-space. 
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Figure  3-17(b)  Ridge  structure  simulation;  formed  by  both  up-dip  and  down-dip 
edges  of  the  basins.  This  ridge  is  constructed  based  upon  the  slopes 
seen  from  the  San  Fernando  and  Los  Angeles  basins.  After  passing 
the  ridge,  surface  waves  are  gradually  developed  from  energy 
trapped  by  the  down-dip  edge  from  Los  Angeles  basin. 
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Angeles  basin  to  produce  surface  waves  at  further  distances.  The  geometry  of 
the  ridge  is  an  important  for  determining  the  amount  of  the  energy  which  is 
trapped  in  the  second  basin. 

Assuming  Heaton’s  source  model  (  1982  ),  the  rupture  delay  between  source 
Si  (  at  13  km  )  and  S2  (  at  7  km  )  is  estimated  to  be  about  4  seconds.  Assuming 
that  Si  and  S3  (  at  1  km  )  represent  areas  of  the  largest  dislocations  ,  we  can 
roughly  estimate  the  timing  of  arrivals  and  corresponding  wave  types  seen  along 
profile  I.  The  time  delay  between  S3  and  SI  due  to  rupture  is  about  6  seconds 
and  the  corresponding  travel  time  delay  is  about  -4  seconds.  Thus  we  concluded 
that  the  S!  phases  observed  along  profile  I  are  mainly  derived  from  Si  (  13  km 
source  ).  The  strong  later  arriving  waves  seen  near  the  Santa  Monica  Mountains 
are  surface  waves  from  both  SI  (  deep  source  )  and  S3  (  shallow  source  ).  For 
the  three  stations  in  the  San  Fernando  Basin.  C04B.  J145  and  HI  15.  the  high 
amplitude  waves  seen  in  the  velocity  waveforms  are  combinations  of  surface 
waves  from  both  deep  and  shallow  sources  and  body  waves  from  the  shallow 
source.  Thus,  the  reason  that  the  S!  phase  maintains  strong  coherence  along 
this  structurally  complex  path  is  that  S!  is  derived  mainly  from  the  deeper  part 
of  the  rupture.  The  S!  travel  path  from  deeper  sources  is  affected  relatively  lit¬ 
tle  by  the  presence  of  local  basins. 

b.  Irregular  Structure  along  Profile  n  Ground  motions  observed  along 
profile  11  show  very  different  features  from  those  observed  along  profile  1.  Again, 
we  generate  synthetic  seismograms  along  the  profile  II  for  three  source  depths, 
namely  1  km  (  source  S3  ),  7km  (  source  S2  ),  and  13  km  (  source  Si  ).  The 
results  are  illustrated  in  Figures  3-18(a),  3-18(b)  and  3-18(c),  respectively. 
Since  the  wave  path  is  more  or  less  following  the  northeastern  boundary  of  the 
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Los  Angeles  Basin,  the  wave  development  is  quite  different  from  that  seen  along 
profile  I.  We  approximate  the  structure  along  this  profile  as  a  basin  (  shear 
velocity  of  2.1  km/sec  )  embedded  in  a  half-sapce  (  3.5  km/sec  ),  as  illustrated 
in  the  left  of  Figure  3*18(a).  In  spite  of  this  rather  simplified  structure  model  , 
the  synthetics  compare  well  with  the  motions  observed  along  profile  11  (  Figure 
3-5). 

As  was  the  case  for  profile  1,  that  the  SI  phase  is  primarily  composed  of 
direct  S-wave  arrivals  from  Si  (  source  at  13  km  ).  In  our  synthetics,  surface 
waves  are  gradually  developed  after  distance  of  about  40  km.  which  agrees  well 
with  the  observations.  According  to  our  2D  modeling,  the  shallow  source  (S3) 
contributes  little  to  ground  motions  observed  along  this  path.  Because  the 
structure  in  the  source  region  is  approximated  by  a  half-space,  little  energy 
from  the  shallow  source  is  trapped  within  the  basin. 

Discussions  and  Conclusions 

Velocity  profiles  recorded  from  the  1971  San  Fernando  earthquake  demon¬ 
strate  strong  path  effects  due  to  topography  of  the  basement  surfaces.  The 
different  waveform  pattern  observed  along  three  long  profiles  implies  that  valley 
or  basin  structures  can  produce  significant  reverberation  and  amplification 
effects.  It  is  suggested  that  the  proper  geometry  of  basins  and  earthquake 
sources  needs  to  be  considered  in  order  to  predict  waveforms  along  these  wave 
paths. 

Velocity  profiles  along  three  local  arrays  in  the  Los  Angles  basin  suggests 
that,  within  3  km  in  range,  the  waveforms  are  almost  identical.  It  is  also  sug¬ 
gested  that  higher  frequency  waves,  as  represented  by  ground  accelerations. 


are  less  affected  by  these  large  scale  basin  structures. 

Profiles  of  displacement  traces,  which  are  not  shown  in  this  study,  display 
characteristics  similar  to  those  observed  in  the  velocity  profiles.  That  is,  the 
waveforms  show  good  correlation  Uh  variations  in  the  basement  surfaces.  How¬ 
ever,  due  to  its  high-frequency  nature,  the  S!  phase  is  more  easily  recognized  in 
velocity  records  than  in  displacement  records. 

Along  profile  III,  displacements  have  relatively  simple  waveform.  This  is 
consistent  with  our  understanding  of  the  basement  surface  along  profile  II,  since 
a  half-space  seems  to  be  an  adequate  approximation  for  the  propagating  path. 
We  speculate  that  only  deeper  sources  effectively  transmit  short-period  energy 
along  this  profile  as  suggested  in  Figures  3-  18(a)  through  3-18(c).  This  is  prob¬ 
ably  the  main  reason  that  the  duration  of  high-frequency  motion  is  rather  short 
for  all  the  stations  recorded  along  this  azimuth. 

We  use  the  2D  finite  difference  method  to  compute  synthetic  seismograms 
for  irregular  structures  which  exist  along  profile  I  and  profile  II.  Although  the 
structure  used  in  the  calculation  is  rather  simplified  and  the  numerical  code  is 
for  acoustical  media,  the  general  features  of  the  synthetics  agree  quite  well  with 
the  observed  ground  motions.  We  approximate  Heaton's  (  1982  )  source  model, 
which  has  two  major  areas  of  dislocation  by  point  sources  ,  at  depths  of  13  km 
and  1  km.  Using  the  relative  timing,  amplitude  and  signed  duration  of  major 
arrivals  observed  in  the  ground  motions  along  both  profile  I  and  profile  II,  we 
suggest  that  the  coherent  S!  phase  (  the  S  waves  )  is  primarily  from  the  first 
source,  at  13  km.  The  shallow  source  (  near  1  km  )  contributes  little  to  the  sta¬ 
tions  along  the  profile  II.  However,  it  does  produce  strong  surface  waves  for  the 
stations  in  the  San  Fernando  Basin  region  and  also  contributes  to  the  generation 
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